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CN ' Abstract 

> : 

("^ , We test the asymptotic scaling behavior of state-of-the-art simulations of QCD with 

two flavors of light Wilson fermions. This is done by matching vr and p masses on 
lattices of size 16^ x 32 and 8^ x 16. We find that at /? = 6/5^ = 5.3 matching 



I is not possible over a range extending down to (3 = 3.5. The large lattice data at 

' P = 5.5 matches the small lattice values at /? = 4.9(1) leading to a shift A/3 = 

0.6(1), considerably larger than the perturbative prediction of 0.45. In both cases 
, we conclude that the simulations are very far from the asymptotic scaling region. 

^ : 1 Introduction 



Assuming that a universal continuum limit of lattice QCD exists, measurements of some 
physical observable Q, such as a hadron mass, in a Monte Carlo computation will not depend 
5^ \ on the details of the lattice discretization if the lattice spacing a is small enough. The value 
of the observable, = a^f{g,K), obtained by numerical calculation in lattice QCD with 
Wilson fermions, depends only on the parameters appearing in the lattice action: the gauge 
coupling (3 = and the hopping parameter k. The dimensionless quantity / has been 
related to the dimensionful observable Q by multiplying by the appropriate power of the 
lattice spacing a. As we approach the continuum limit a — the fact that Q is independent 
of a is expressed by the renormalization group equation 

og OK, J 

where 

n dg dn 
p = —a—, p = —a—, 
da da 

The derivatives with respect to a are taken at "constant physics." As we have two in- 
dependent relevant parameters we require two independent physical observables to use as 
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renormalization conditions to fix these parameters; j3 and (3' should, of course, be indepen- 
dent of the choice of renormahzation conditions. 

In a perturbative expansion of (3 about (7 = only the first two terms are independent of 
the details of the lattice regularization scheme. When the higher order terms are negligible 
one generally speaks of being in the asymptotic scaling regime. The aim of this work is to 
determine, for the realistic case of QCD with two flavors of light Wilson fermions, how close 
present simulations are to this asymptotic scaling regime. 

In numerical simulations it is impractical to determine the /5-function, which correspond 
to an infinitesimal change of scale, directly. Rather we work with a finite change of scale. 
We consider two lattices, Ai with Li sites in each direction and with lattice spacing oi, and 
A2 with L2 sites in each direction and lattice spacing a2- By fiat we assert that Ai and A2 
represent the same physical volume, thus LiQi = L2a2- 

If we select some set of physical observables fi, we have 

where 5i is the dimension of Vti and fi is the corresponding dimensionless function measured 
on the lattice. In the continuum limit all these observables must become independent of the 
lattice, so for the two different lattices 

We define the quantities 

. „ 6 6 . / \ / \ 

= -IT, T TT, T, A/t = K[a2) - K[ai) 

to be the change in the couplings needed to compensate for this change in the cutoff. 



2 The /3-function 

We expect the /9-function to become independent of k in the continuum and chiral limits.0 
The /3- function at fixed k is defined by (3{g) = —adg/da, and the two universal terms 
from perturbation theory are f3{g) = —h^g^ — hig^ + 0{g'^) where 



33 - 2n 



f 



487r2 
153 - 19nf 

3847r4 



bo = 
bi = 

for QCD with n j flavors of light fermions. 

^If we use continuum perturbation theory with dimensional regularization and the minimal subtraction 
renormalization scheme, then the BRS identities insure that all divergences are logarithmic, and thus must 
be independent of the quark mass. From this we may conclude that in the vicinity of the chirally symmetric 
continuum limit any k dependence of the /3-function must arise from irrelevant lattice operators and be 
suppressed by some power of the lattice spacing. 
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The solution of these equations, expressed in terms of /3 = 6/(7(02)^, is 

In _ 76 /"^ dx 

after evaluating the integral we find that 

A/5p,rturbativc ^-12 ^60 + ^^ In (^^^ 

which is the perturbative prediction for the change in j3 needed to effect a change in scale 
by a factor 02/01. 

A/9 has been thoroughly studied for pure gauge theory with the Wilson action (see [Q 
and references therein) and exhibits the well-known dip around (3 = 6.0. For QCD with 
dynamical fermions we are unaware of recent studies of A/5. For staggered fermions, Blum 
et al. [0, determined the /5-function, needed for the computation of the equation of state 
for finite temperature QCD, from fits to hadron masses obtained at different couplings, and 
hence different lattice scales. They found significant deviations from asymptotic scaling that 
again would lead to a dip in A/?. 



3 Results 

The /5-function for two flavors of Wilson fermions was computed earlier by some of us 0, 
The physical observables used for that purpose were m^r and several Creutz ratios of Wilson 
loops. 

The main conclusion then was that perturbative scaling was only expected to occur for 
f3 > 6.0. In those simulations the vr-like object was heavy, and it was thought at the time 
that this may have hidden the effects of the dynamical fermions in the system. Simulations 
with smaller 771^^ are thus needed to clarify the issue. 

In this work we concentrated on data available on lattices of size 16^ x 32 with two flavors 
of Wilson fermions with parameters and results shown in Table ^ 0. These were matched 
using 8^ X 16 lattices with the parameters and results for the interesting parameter region 
shown in Table |^. 

A preliminary report on this work was presented in 0. 

The values of f3 were chosen for hadron spectrum and other calculations largely based 
on the assumption that the presence of two fermion flavors renormalizes /3 downwards by 
about 0.5 from its quenched value, and thus we may not be far from the perturbative scaling 
regime at these parameter values. 

For our test of scaling we used two hadron masses, ttt-tt and nip, as observables for matching 
with similar data generated on lattices of size 8^ x 16 at values of (3 ranging from 5.1 down 
to 3.5 and many k values. 

The 8^ simulations used step sizes from 0.015 to 0.025 to insure a reasonable acceptance 
rate. The conjugate gradient residual used for the molecular dynamics steps was 10~^. 
Hadron measurements were made every 5 trajectories. We found integrated autocorrelation 
times were typically of this order. 
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/5 




# 


A 




arup 


arriN 


m^/nip 




5.3 


0.1670 


484 


5 


0.4540(20) 


0.6350(20) 


0.9620(40) 


0.715(5) 


0.4719(29) 


5.3 


0.1675 


417 


3 


0.3120(40) 


0.5230(40) 


0.7660(90) 


0.596(15) 


0.4073(71) 


5.5 


0.1596 


400 


5 


0.3754(27) 


0.4947(78) 


0.7644(48) 


0.759(11) 


0.4911(36) 


5.5 


0.1600 


400 


5 


0.3262(17) 


0.4606(28) 


0.7106(91) 


0.708(4) 


0.4591(57) 


5.5 


0.1604 


669 


3 


0.2666(20) 


0.4208(55) 


0.6326(56) 


0.634(8) 


0.4215(48) 



Table 1: Parameters and masses on 16'^ x 32 lattice [^,0. # is the number of measurements 
made, each separated by A trajectories. 







# 




arrip 




m^/mp 




4.8 


0.1^ 


^90 


130 


0.6681(41) 


0.964(12) 


0.964(12) 


0.678(16) 


0.4038(51) 


4.8 


0.1900 


170 


0.5814(35) 


0.897(7) 


1.522(26) 


0.6480(54) 


0.3820(64) 


4.8 


0.1905 


109 


0.5110(48) 


0.896(28) 


1.362(26) 


0.571(19) 


0.3753(65) 


4.9 


0.1^ 


m 


189 


0.7826(28) 


1.013(5) 


1.783(36) 


0.7722(36) 


0.4390(88) 


4.9 


0.1^ 


m 


190 


0.7520(23) 


0.988(5) 


1.714(20) 


0.7609(33) 


0.4388(49) 


4.9 


0.1^ 


^50 


110 


0.6808(46) 


0.928(5) 


1.549(13) 


0.7334(44) 


0.4396(40) 


4.9 


0.1^ 


^55 


272 


0.6491(40) 


0.930(7) 


1.660(57) 


0.6982(63) 


0.391(13) 


4.9 


0.1^ 


^60 


200 


0.5636(47) 


0.872(10) 


1.451(31) 


0.6467(88) 


0.3886(88) 


4.9 


0.1^ 


^65 


150 


0.5018(59) 


0.852(20) 


1.55(11) 


0.590(13) 


0.324(22) 


4.9 


0.1^ 


^70 


70 


0.4981(65) 


0.816(20) 


1.395(23) 


0.611(15) 


0.3572(81) 


4.95 


0.1^ 


^15 


80 


0.8083(37) 


1.003(6) 


1.693(19) 


0.8059(44) 


0.4775(47) 


5.0 


0.1^ 


^00 


30 


0.7734(85) 


0.990(14) 


1.690(36) 


0.781(11) 


0.4578(95) 


5.0 


0.1^ 


^10 


150 


0.6321(53) 


0.875(7) 


1.462(19) 


0.7228(53) 


0.4325(65) 


5.0 


0.1^ 


lib 


270 


0.5145(83) 


0.814(9) 


1.342(19) 


0.632(12) 


0.3834(77) 



Table 2: Parameters and masses for the matching region on 8^ x 16 lattices. 7^ is the number 
of measurements made. 
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Quark propagators with wall sources and point sinks were computed in Coulomb gauge 
on time slices and Lt/2. We determined the meson masses by a correlated single state 
cosh fit to two correlation functions with interpolating fields il)Tiip and '?/'7orj'?/;; our baryon 
correlation functions used only one interpolating field. The covariance and errors were de- 
termined via the bootstrap method. Fits were made through the center of the lattice with 
a starting time slice determined by the criterion max A^^ofQ 5m/ where N^^f is the number 
of degrees of freedom for a fit of mass m. In all cases Q, the 1 standard deviation confidence 
level, was greater than 10%. 



3.1 Non-Matching at = 5.3 

Our first result is the unmatchability of the data at /3 = 5.3. We find that A/? must be larger 
than 1.8. Figure |l| shows the two points obtained at /? = 5.3 on the larger lattice (with m.„ 
appropriately scaled by a factor of 2) and the corresponding data for and m-^/mp on the 
smaller volume for P = 3.5. Data points obtained at larger values lie to the left of those 
from P = 3.5 (see figure 1). 



0.90 




Figure 1: Non-matching of 771^^/ nip at /? = 5.3. 
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3.2 Matching at /3 = 5.5 

Our second result is that at /3 = 5.5 matching is possible. Figure |^ shows the matching with 
(3 = 4.9 on the smaller lattice. It is interesting to note that all three points, at different 
values of n, match with corresponding points at the same value of (3. For comparison similar 
points obtained at jS = 4.8 and 5.0 on the smaller lattice are shown. For (3 = 4.8 mT^/rrip 
tends to be lower than the values at /? = 5.5, whereas at = 5.0 they are higher. Two 
conclusions follow from this data: the first is that the data at /3 = 5.5 is consistent with a 
value of A(3 = 0.6(1), and second that this matching is the same for all three different values 
of K used, indicating a k independence of the /3-function. 
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Figure 2: Matching of m^r/mp at (3 = 5.5 

At the central matching (3 value {(3 = 4.90) we find Ak = -0.0249(1), -0.0253(2), and 
—0.0258(2) for large lattice values k = 0.1596, 0.1600, and 0.1604 respectively, thus giving 
an estimate of the integrated /^'-function. 

On the other hand, figure |^ shows that the the system is not in the scaling region because 
the pion-nucleon mass ratio does not match for the same values of (3 and k. 

Another set of observables which we can match to determine Af3 are appropriate ratios 
of Wilson loops, just as was done for pure gauge theories p. We computed space-like planar 
Wilson loops on all our lattice configurations and evaluated ratios R of two or four Wilson 
loops arranged such that the corner and perimeter singularities cancel between numerator 
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Figure 3: Non-matching of uiT^/mN aX (3 = 5.5 

and denominator, 

R{l3,K,{ri,r2,r3,r4},Ls,Lt) = — r, n + r2 = + r^, 
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(ra, r4) 
W{ri,r2)W{r3,r^) 
W{r5,r6)W{r7,r8) 



i=l i=5 



We then determined A/5 from the condition 

RiP - A/3, k', {r}, iL„ iLi) = i?(/?, k, {2r}, L„ L*). 

on the smaller lattice was tuned such that the vr mass on the smaller lattice was twice 
that on the larger lattice. In practice we interpolated (or extrapolated) linearly from two 
simulations on the smaller lattice with slightly different values of (3'. Unfortunately, with 
our statistics the errors on the larger Wilson loops grew quickly and they were not useful for 
the ratio matching. With smaller Wilson loops dominating the ratios lattice artefacts in the 
matching become important. We tried to avoid this by considering tree level and one- loop 
improved ratios [§]. For the matching of the k = 0.16 data, for example, depending on various 
"cuts" imposed on the ratios considered, such as the minimum perimeter and area of the 
Wilson loops and the maximum relative error, we obtained average AjS values between 0.36 
and 0.55. This is somewhat smaller than the preferred value from the matching of m.,t/mp, 
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indicating again that we are not yet in the regime of a unique, observable independent 
/^-function. 

The static quark potential provides us with another matching condition between the 
large and small lattices. We computed finite t approximations to the static quark potential 
using time-like Wilson loops W{f,t) which were constructed using "APE"-smeared spatial 
links IP, ^. On- and off-axis spatial paths were used with distances r = n, \/2n, \/3n, and 
\^n, with n a positive integer. The "effective" potentials 



W{f,t+1) 

were fitted using a correlated procedure with the spatial covariance estimated using the 
bootstrap method. The potential on the large lattice was fitted to the formQ 

Vir) = Vo + ar--- e'(Gi(f) - 

the last term takes account of the lattice artifacts at short distance. Here, Giir) denotes 
the lattice Coulomb potential for the Wilson gluonic propagator. 

In figure ^, we plot the potential for the 16^ x 32 lattice with j3 = 5.5 and k, = 0.1600 
scaled to a spatial lattice of length 8. Also plotted are the potentials for the 8'^ x 16 lattices 
with f3 = 4.80 and k = 0.1890, and (3 = 5.0 and k = 0.1810: the latter point only has the 
on-axis values shown. These points have nearly matching pion-rho ratios. 

We can see that the 8'^ x 16 potentials do not match the 16^ x 32 lattice. While we do 
not expect the constant Vq to match between the fine and coarse lattices, it is clear that 
the slopes of the coarse lattice potentials are too large. The slope decreases with increasing 
/3, indicating that the shift A/5 required for matching the potential is smaller than that 
required for matching m^/mp and closer to that of the Wilson loops ratios. We also see 
large violations of rotational symmetry in the coarse lattice potential compared to the finer 
lattice potential, indicating that the coarse lattice is at relatively strong coupling. 

3.3 Parameter reduction at strong coupling 

A third result is an approximate parameter reduction in the theory at strong coupling (/3 = 
3.5 to 13 = 4.7). Figure ^ includes all our data in this range of /3 on the smaller lattice and 
shows that for 3.5 ^ (3 < 4.9 the data is essentially a function of only a single parameter. 
Any change in f3 may be compensated for by a change in k. This is consistent with the large 
N strong coupling result of Kawamoto and Smit [|l^ , which for = 3 and r = 1 gives 



ml = 4.8m, = 0.984 + 1.97m. 



4 Conclusions 

Combining the present results with those obtained earlier on smaller lattices we find 
the /3-function shown in figure |[ It is similar to that found for pure SU{3) lattice gauge 

^We note that we never observed string breaking. 
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Figure 4: Non-matching of the qq potential The hne shows the "effective" potential fitted to 
the 16^ X 32 lattice potential at /3 = 5.5, with the solid part of the line indicating the range 
of r values used in the fitting procedure. The best fit gives Vq = 0.748(4), a = 0.0585(6), 
e = 0.340(6), and e' = 0.33(1) with = 8-7 for 14 degrees of freedom. 

theory with the Wilson gauge action. We are led to the conclusion that current dynamical 
simulations are being done at strong coupling. 

Indeed, not only do we find that for j3 < 5.5 the /5-function does not agree with the 
perturbative predictions (asymptotic scaling), but that different renormalization conditions 
give different estimates for the /9-function: the shift Af3 = 0.6(1) which gives satisfactory 
matching for m^r and iriT^/mp fails to match m.,^/mN, ratios of Wilson loops, or the qq 
potential. In other words we do not even observe a universal but non-perturbative /3-function 
(scaling) . 

One therefore may need to reassess estimates for parameters and hence the cost of future 
computations with dynamical Wilson fermions accordingly. 
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Figure 5: Cumulative data for 3.5 < (3 < 5.1 on an 8^ x 16 lattice. The line shows the large 
strong coupling result. 
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